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ABSTRACT 

The surface of a large, anechoic tank was excited-to produce several 
different anisotropic sea and swell surfaces characterized-by “typical” 
but distinctive frequency spectra and wave height probability densitiés. 
Surface roughness for specular scatter, g (=4k“o%sin“e), was varied 
from zero to four by selection of underwater sound propagation constant, 
k, RMS wave height, o, and grazing angle, 6. Concurrent measurement 
was made of 20 second averages of squared surface wave height, and 
coherent specularly-scattered sound intensity. For g < 1., coherent 
specular scatter from all surfaces, including pure swell (with a 
strongly non-Gaussian wave height probability density) was found to vary 
as e 9, as previously predicted for Gaussian surfaces only. For certain 
combinations of 98, co, and sound wavelength, 21, coherent components 
significantly larger than predicted were noted in both cross-wind 
and down-wind scattering for g > 1.3; this is believed to be a previously- 
hypothesized interference phenomenon. Cross-sea scattering consistently 
produced a slightly stronger coherent component of intensity .than 


scattering in the direction of the sea. 
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I. INTRODUCTION 


The last ten years have witnessed the development and deployment of 
low-frequency, high-power, surface ship sonars. Intelligent use of 
these long-range sensors, and the ability to predict their performance 
under a given set of conditions, require quantitative understanding of 
the sound propagation involved. Both surface and seabed reflections of 
the beam can greatly affect sonar performance; quantitative analysis of 
“bottom-bounce" phenomena generally has not proved fruitful because of 
the dearth of knowledge concerning the geography, topography, composi- 
tion, and acoustic properties of the ocean bottom. Analysis of surface 
scattering, however, is more tractable because of the comparative ease 
with which the properties of the reflecting boundary may be studied. 

In addition, surface scattering experiments are easily scaled to labora- 
tory dimensions. 

When a beam of underwater sound is incident upon the water's surface 
it is scattered, the nature of this scattering being dependent upon the 
characteristics of the interface. If the surface is smooth, the scat- 
tering will be purely specular. Conversely, if the surface is rough, 
sound will be scattered in all directions, with signal strength being a 
function of the angle of observation. Rough-surface reflection produces 
a scattered sound field which contains both coherent and incoherent pres- 
Sure components relative to the incident pressure. The behavior of the 
coherent, scattered component is of prime importance because of the in- 
creasingly more common use cf correlation detectors in sonar systems. 

Scheible and Fowler (1968) studied snecular scatter of underwater 


sound from a wind driven model sea at the Naval Postgraduate School; 


their investigation included both the amplitude distribution and the 
frequency spectrum of the scattered sound. 

The work reported herein is essentially an extension and refinement 
of the work done by Scheible. An analogue correlator (which Scheible 
did not have) was used to isolate the coherent component of the scat- 
tered sound pressure. As will be explained in a later section, this 
isolation permitted the use of much "rougher" water surfaces for scat- 
tering. The types of scattering surfaces used were increased to include 
pure swell, sea-swell combination, and wind driven seas of three differ- 
ent heights. In addition, scattering both parallel and perpendicular to 
the sea direction was studied. 

For all of these "seas" the objective has been to use the model ex- 
periment to study cases that are beyond the capability of theorists, and 


beyond the practicality of ocean-acoustic experiments. 


II. THEORETICAL BEHAVIOR OF THE SCATTERED SOUND 


A. THE SPECULARLY SCATTERED INTENSITY 

The following treatment briefly summarizes the statistical approach 
of Medwin [Ref. 1], which in turn follows the development of Beckmann 
[Ref. 2]. The measure of the roughness of the water surface from which 


the sound is scattered is defined by the roughness parameter g. 


be oO Z 
g = [2n (cos o, + cos 5) ] 


where: o = Root-mean-squared (RMS) Wave Height 
A} = Wavelength of sound in water 
o, = Angle of incidence of sound upon the surface 


(measured with normal to surface) 


bo = Angle of reflection of the sound 


For specular scatter, o, = bo = 4 and g becomes 


g = k%o“(2 cos 4)° where: k = 2 = > 
= 4h¢+careyp 
The general scattering problem is then attacked by applying the 
Helmholtz Integral to a volume within the water, part of the surface of 
which is the ensonified surface area. If the Kirchhoff assumption (the 
use of the plane wave reflection coefficient at each point on the sur- 
face) is satisfied, an expression for the time-averaged, reflected 
intensity relative to mirror surface intensity results, which, for the 
Specular direction reduces to: 
I 
ee = e@ 9 + £(A,L,¢,0,A) 


mirror 
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where L = Correlation distance on the rough surface 


I> 
I 


Ensonified surface area 


I = p*p where p = acoustic pressure 


The above expression was derived for a surface correlation function of 
Gaussian form, and for a Gaussian distribution of surface wave heights. 

The first term, describing the coherent specular component, is iden- 
tical to that obtained by several authors. The second term represents 
Scattered sound which is incoherent with the transmitted signal; since 
it was not possible to experimentally isolate the very small incoherent 
component, the exact form of the term is not given here. However, the 
following theoretical conclusions are considered worthy of note: 

1. The magnitude of the incoherent term is a function of ensonified 
area and correlation distance, while the coherent term, eS, 1s 
incenence of these parameters. 

2. In the mirror surface limit (g>0) the incoherent term vanishes, 
and the coherent component approaches unity. 

Scheible [Ref. 3] used this limiting behavior to suppress the incoherent 
term; he limited g to 0 < g < 0.3, stated that the first term was much 
greater than the second, and dropped the incoherent term. 

Attacking the problem through the surface wave spectrum, Parkins 
[Ref. 4] also predicts a coherent term of the form e 9. His incoherent 
term arises from selective interactions between appropriate surface wave 
and sound wavelengths. In the purely specular direction, zero incoherent 
component is predicted, since his treatment considers only the central ray 
and neglects scattering contributions from off-axis parts of the ensonified 


ared. 


B. SURFACE SHADOWING AT NEAR GRAZING INCIDENCE 

As the grazing angle (complement of angle of incidence, 9) approaches 
zero, shadowing of certain areas of the scattering surface will result. 
A shadowing correction, derived from ray acoustics and a statistical 
representation of the surface, is presented by Wagner [Ref. 5], who 
states that the correct use of the shadowing function, S(6), for the 


coherent scattered component is: 


I 
coherent : 52(9) ae 
mirror 
where: S(e) = (1 - e~ "By erf (zy) 
-1/2 _-(9v*)/8 -1/2 _-y" 
B = x L(3nv) e + (nv) ““" e - erfc v] 
—. tan 6 
2(o/L) 
o = RMS wave height 
L = Surface correlation distance 
6 = Grazing angle 


Diffraction effects are neglected by the ray treatment; it is therefore 
expected that the accuracy of Wagner's correction will increase with 


increasing aan 


III. EXPERIMENTAL PROCEDURE 


A. THE ANECHOIC TANK 

The scattering experiments were carried out in the north anechoic 
tank, filled with fresh water, and of dimensions 24' long, 6' wide, and 
8' deep. A wall and bottom lining of acoustically absorbant rubber 
material renders the tank essentially anechoic at frequencies above 
20 kHz. A swell maker and wind-wave generator were located at the east 
end of the tank as shown in figure 1. Clear plastic covers, suspended 
6" above the water, confined the wind generator's output to the surface, 
producing a wind-driven sea. Five separately controllable, centrifugal 
fans feeding a common plenum provided the flexibility needed to nnoewee 
RMS wave hevante non 0.1 to 0.5 cm. Swell waves were produced by sinu- 
soidally oscillating a "V" shaped wooden wedge in the vertical plane by 
means of a motor-driven eccentric, as shown in figure 2. 

The sea, swell, or sea-swell combination propagated down the tank, 
past the sound scattering area, and was largely dissipated by a sloping 
"beach" consisting of loosely packed aluminum shavings within a plastic 
screen enclosure. In all cases the ensonified area was centered on the 
wave staff, which was located at a fetch distance of 3.32 meters meas- 
ured from the plenum discharge. 

The surface area ensonified by the active transducer can be computed 
geometrically by assuming no ensonification outside the -3dB angles. The 
resulting surface is an ellipse, whose area, A, and largest linear dimen- 
sion, d, are computed as shown in figure 3. Attention is invited to the 
distinction between grazing angle, 6, and its complement, angle of inci- 
dence, ¢. The grazing angle, 6, was determined by measuring transducer 


separation and transducer depth. 
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To scatter sound both parallel and perpendicular to the sea direc- 
tion, three separate transducer arrangements were used as shown in 
figure 4. The transducers, supported by the noeMeinawracks described 
in Ref. 1, were always positioned with their axes in the specular direc- 
tion. For frequencies up to and including 120 kHz, two EDO model 327 
piezoelectric transducers were employed; the face of each consisted of 
two concentric, electrically separated driving sections. At a given 
frequency, two beam patterns were available: a broader beam ob tanned 
by driving the inner section only, and a narrower beam resulting from 


driving both sections in electrical parallel. 


B. THE SEA SURFACE 
Various surfaces consisting of "sea" alone, "swell" alone, and 
“sea-swell'" combination were used in the scattering study. To determine 
the nature of each scattering surface the instantaneous wave height, 
z(t), at the center of the ensonified area was continuously monitored by 
an omnidirectional wave staff; the resulting wave record was tape recorded. 
Upon completion of a particular scattering study, the corresponding wave 
data were replayed and processed to determine the RMS wave height which 
existed at the time of the sound scattering. 
1. The Wave Staff 
The wave staff and its associated electrical circuitry was de- 
Signed, fabricated, and calibrated locally. The staff, shown in figure 5, 
is of the resistance type and operates as one leg of a balanced resistance 
bridge, as shown in figure 6. 
Over a limited range of immersion the resistance between the staff 
and ground rod is inversely proportional to the depth of immersion, and it 
is within this linear range that the staff is made to operate, as shown in 


figure 9A. 
18 
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Preliminary experimentation revealed that the position of the 
around rod relative to the staff affected the static calibration of the 
instrument. To eliminate this variable, and to insure its insensitivity 
to passing waves, the ground rod was mounted on the same support bracket 
as the staff and was insulated by a plastic sleeve where it pierced the 
water's surface (figure 5). 

The wave staff assembly formed a time-varying resistance whose mean 
value was about 2500 ohms; the other three legs of the nearly-balanced re- 
sistance bridge, shown in figure 6, were 2500 ohm resistors. This bridge 
was fed with AC voltage at 25 kHz to prevent the formation of electrolysis 
bubbles on the staff, which would have caused a change in resistance. The 
bridge output, a 25 kHz voltage with a time-varying envelope, was amplified, 
demodulated, and reduced to zero mean by a high-pass RC filter with half- 
power point of 0.05 Hz. The output of the circuit is called the wave record, 
z(t), a time-varying voltage proportional to the instantaneous wave height. 

a. Wave Staff Calibration 

The data of figure 9A (from which was computed the static probe 
calibration) were obtained by measuring changes in probe circuit output 
voltage as the water level in the tank was changed by accurately measured 
increments; it was necessary to take these readings with the DC blocking 
capacitor of figure 6 by-passed. 

In addition to this static deflection calibration, the wave 
staff (and its associated circuitry) was calibrated for surface wave fre- 
quencies up to 30.0 Hz. At each frequency checked, an independent 
determination of the calibration was made. Each was then normalized 
(figure 9B) by dividing by the previously obtained static calibration. 

When the resulting data were plotted as shown in figure 9B, the frequency 


response curve for the probe and its circuitry results. In the frequency 
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range 0.5 < f < 8.0 Hz, motion of the water relative to the staff was 
obtained by vertically oscillating the staff assembly in a beaker of. 
"tank water", as shown in figure 7. The frequency and amplitude of this 
displacement were both adjustable and accurately measurable. 

In the frequency range 8.0 < f < 30.0 Hz, change in probe 
immersion was accomplished by generating capillary waves on the surface 
of a small] tank. Referrring to figure 8, the generating mechanism was 
a lightweight fiber cone, mounted point down, partially beneath the sur- 
face as shown. The cone was driven sinusoidally in the vertical plane 
by an electric brated (a permanent magnet, moving coil device) mounted 
rigidly above the tank. The frequency and amplitude of the cone's 
motion could be varied by changing the corresponding properties of the 
driving voltage. To prevent capillary reflection from the walls of the 
| rectangular tank, a 5 onal "beach" of blotting paper was installed 
along all four sides. Capillary height was measured both by the wave 
Staff ard by a micrometer assembly, mounted rigidly to the frame of the 
tank, by means of which a thin, pointed rod could be lowered to touch the 
water's surface. Measurement of maximum wave height by the micrometer 
Pee erbily was done as follows: 

1. With the surface in a mirror condition, the rod was 
slowly and carefully lowered until its sical touched the surface. 
Reading the micrometer dial gave a surface height reading eefaté to 
t.02 mm; to reduce this uncertainty, five consecutive readings were taken 
and averaged. 

2. Capillary waves of the desired frequency were then gener- 
ated, and the above measurement repeated. The point of first contact in 


this case was the height of the crests of the waves. 
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3. Capillary height (crest-to-trough) was then computed as 
being twice the difference in the two measurements. This step assumes 
that crest-to-trough height is equal to twice the value of mirror-to- 
crest height. 

The crest-to-trough voltage variation from the wave staff was 
divided into the measured wave height to obtain the staff calibration at 
that particular frequency. This calibration was then normalized by dividing 
by the static calibration. Figure 9 summarizes the staff calibration from 
zero to 30 Hz. The error flags become larger at higher frequencies, be- 
cause wave height decreases rapidly with increasing f, while the absolute 
possible error is a constant of the micrometer scale. 

2. Processing the Wave Record 
a. Temporal Variation of Mean Squared Wave Height 

The correlator, as wired and programmed in figure 10, accepted 
the wave record, z(t), squared it, and averaged the resulting product for 
20 seconds on a continuous basis. The correlator output was a voltage 


representing mean squared wave height, a 


ney = a ee () 
This voltage, displayed on a strip chart, revealed that oa itself was 
itself a function of time, e<(t), in spite of the 20 second integration; 
variations on the order of 20% were observed over a ten minute interval 
and were present for all wind-driven surfaces. The nature of these low- 
frequency temporal variations was not known; they appeared to be statis- 
tical, but on the possibility that there might-have been an external 


cause, several attempts were made to study this effect. 
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Close examination of the strip chart revealed that the fluc- 
tuations in eat were of very low frequency, apparently less than 0.1 
Hz (this observation being consistent with the 20 second time constant 
of the correlator). The possibility that the fluctuations were being 
introduced in the electronics was checked by oscillating the wave staff 
in a beaker of tank water as shown in figure 7. The amplitude of the 
oscillatory motion was a known constant; the resulting function, a when 
displayed on a strip chart was also constant, indicating that the previ-. 
ously observed variations were a real, if unexplained, phenomenon of the 
wind driven surfaces. The electrical power drawn by each of the five 
centrifugal fans was monitored and found to be constant with time. 
As a second possibility, the period of the tank's natural 

‘seiche was determined by shocking the system (quickly opening and closing 
the large drain valve at the tank's west end). The period of the resulting 
water surface oscillation was 3.7 seconds — at least one order of magnitude 
smaller than the periods of fluctuation observed in o@(t). | 

The energy spectral density of the fluctuations was determined 
using the equipment of figure 11. To shift the low frequencies into a 
range more tractable to analysis, the function o“(t) was tape recorded at 
speed 0.375 IPS, and replayed for processing at 37.5 IPS. It was found 
necessary to eliminate the DC component of the signal, sires this would 
produce delta function-like behavior at f=0 on the spectrum plot, possibly 
masing the low frequencies of interest. DC elimination was done with a 
high-pass, RC filter, whose parameters are shown in figure 11. The survi- 
ving signal after this filter was o“ (100t*) - o“(100t*) » (where t* = 
time in the playback frame of reference). The autocorrelation of this 


function was computed by the PAR-101, then the cosine transform of the 
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autocorrelation was determined by the PAR-102, which yielded the energy 
spectral density in accordance with the Weiner-Kinchine Theorem. | : 

The resulting plot, shown in figure 11, revealed three strong, 
low-frequency components with periods of approximately 500, 100 and 50 
seconds. It is possible that these components might have been introduced 
by the high-pass filter, but it is believed that the fluctuations are 
statistical, possibly introduced by non-linear surface wave interaction. 
Because of the behavior of this variable, it was dated to monitor it 
continuously during all scattering experiments, and to record its value 
at the instant that sound scattering data were taken. This procedure. 
appeared both reasonable and consistent, because both 57 (t) and the scat-— 
ceren sound data were processed through integrators with identical time 
constants. THe mechanics of data taking are described in section III-D. 

b. Water Surfaces Used for Scattering 

Five water surfaces were chosen from which to scatter sound, 
covering the widest possible range of RMS wave height, o, and containing 
both swell and sea=swell combination; a tabular summary of these surfaces 
is given in figure 12. Using the sais ete Fourier analyzer, and 
plotter (but not the filter) of figure 11, energy spectral density plots 
were made for the surfaces and are presented in figures 13 through 15. 
For consecutive runs under identical wind source conditions, slight 
shifting of peak size and location was always noted, as well as changes 
in the total area under the curve, a fact previously established since 
area eee the curve is equal to mean squared wave height, o“(t). Figures 
13 through 15 are representitive energy spectra, each compiled from exam- 
ination of several spectra for nominally the "same" statistical surface. 


It should be noted that nonlinear interaction caused the sea and swell 
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Spectrum "E" to be significantly different from superposition of the two 
constituent spectral (sea "B" and swell "D"). 

A closer analysis of the energy spectral density of wind 
driven sea "B" was made by plotting the data on log-log axes. A continu- 
ous, four hour wave record, z(t), was tape recorded, and replayed at 
one-hundred-times speed, 2(100t*) , into a General Radio Wave Analyzer, 
model 1900-A, using a bandwidth of 3 Hz; the data obtained are plotted in 
figure 16. The curve labelled "noise" represents total system noise, and 
was obtained by frequency analyzing a four hour tape of the "wave record" 
of a mirror surface, in the manner described above. The dependence of 
energy upon frequency at the high-frequency dropoff was found to be care 
which is in close agreement with the Neumann-Pierson spectrum for wind 
waves and Phillips theory [Ref. 6] which predicts a Scheible [Ref. 3] 


obtained fr2+9 


for a slightly smaller sea in the same tank. 
c. Inhomogeneity of the Water Surface 

The wind driven seas produced were not statistically homogene- 
ous. The degree of non-homogeneity in the RMS wave height was determined 
by analysis of 30 minute wave records from the nine surface locations 
shown in figure 17. A wave record from any one of the nine locations was 
processed as shown in figure 10. The output of the correlator was a 30- 
minute record of ote which was tape recorded at .375 IPS. This record 
was then replayed at 3.75 IPS into the analogue-tc~-digital converter of 
the SDS 9300 digital computer. The converter sampled the input function, 
o“(10t*), at a 30 Hz rate. The square root of each sample was then ex- 
tracted, which resulted in an ensemble of approximately 5400 digital 
values of the RMS wave height, o. The mean and standard deviation were 


computed for each of the nine ensembles, and the results are presented in 
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FIGURE 17 INHOMOGENEITY OF THE 
RMS WAVE. HEIGHT OF WIND 
DRIVEN SURFACE "B” 
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figure 17. Wind driven sea "B" was used throughout; a fetch-limited 
Situation is indicated by the significant increase of wave height with 
fetch distance. 
d. Wave Height Distributions for the Scattering Surfaces 

As previously mentioned in section II-A, the e~9 term for 
coherent specular scatter was derived by Beckmann (and all others) after 
assuming a Gaussian distribution of wave heights [Ref. 2]. Wave height 
distributions were determined for the various scattering surfaces and 
are presented in Appendix B, along with a detailed description of how 
the data were obtained. In summary of this appendix, it can be stated 
that all wind driven seas — including the sea-swell combination — had 
height distributions which resembled the Gaussian, although all were 
skewed to the right. The pure swell, however, had a height distribution 


which resembled that of a sine wave. 


C. ANALYSIS OF THE SCATTERED SOUND PRESSURE 

Continuous-wave, rather than pulsed, sound was used in the scattering 
study; however, at each frequency and for each transducer position, 
pulsed sound was used to determine the magnitude of the unwanted direct 
ray between sender and receiver, as well as other non-surface reflections. 
In most cases the ratio of direct=-path to mirror-reflected pressure was 
about .01, and in no case did it exceed .03. Unwanted reflections from 
the tank walls were generally one order of magnitude smaller than the 


direct path signal. 


1. Isolating the Coherent, Received Sound Pressure 


’ 


It was assumed that the sound pressure, p, received at the hydro- 
phone after rough surface reflection was composed of two components, one 


coherent with the transmitted signal and the other incoherent. 
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a Pooh e Pincoh 


Use of the PAR Correlation Function Computer made possible the isolation 
of the coherent component, and (in principle) determination of the magni- 
tude of the incoherent component. 

The acoustic pressure corresponding to the driving voltage was 
represented by pyft) = PA cos(wt). When the surface was in the mirror 
condition, the output of the receiving hydrophone was totally coherent 
with the sent signal. After amplification this received signal was repre- 
sented by BP cos w(t-t'), where B is an amplitude scale factor, and wt' 
represents the phase shift due to propagation distance and possible phase 
Shifts within the transducers. Cross-correlation of sent and received 
Signals produced the function M(t). 

B 20 sec 

M(t) = —=7y ; if Py cos w(t-t ) Cos a t=t*) dt 
where t = correlation delay applied to the pylt) input. 

The output of the correlator was displayed on an oscilloscope, from which 
the value of M(t) could be read as a function of the variable t. The 
function has a maximum at t= t'. 


B 20 sec 


_ ; Bs Z 2 , 
Meee = M(t') = ay ‘ f Py cos w(t-t') dt 


2 


= B Py cos“(wt-wt") = is py (t) =@B . I 


d 


Ly iS proportional to driving intensity, and () means an average over 
20 seconds. | 

The water surface was then roughened and the correlation repeated, 
taking care not to change the driving pressure, pyit)- The scattered 


pressure was now of the form 


4| 


p = ABP,cos (t-t') + p(t) 


d incoh 


where A = amplitude reduction factor due to rough surface. Cross- 
correlation of sent and received signals then gave N(t), which was composed 


of two integrals. 


AB 20 sec 0 

N(t) = 30 : f Py cos w(t-t') cos w(t-t) dt in 
Py 20 sec 
a) ; f BX Caney - COS w(t-t) dt 


The second term vanishes because of the incoherence of the two factors 


in the integrand. The first term has a maximum at t= t'. 


AB 20 sec 9 


z ty 8 2 ; 
N cae =e (dene = 50 ; f Py cos”) w(t-t') dt 
: 2 - 
= AB Py (t) = AB > 
Forming the ratio of Ne to Pl, gives: 
ce 
N(x) aa de a 
= te) 
B /I 
d 
Since A relates pressures, the intensities are related by Ac: 
I(coh) 2 
I(mirror) 


4 - 2 
where (eon) 7 (PconP*con and Lnirror = 8B ly 


The theory of section II-A predicts that ae =e 9, 


2. The Incoherent Component of Scattered Sound Pressure 


In principle at least, the ratio of incoherent to mirror intensity 


42 





can be determined using the correlator in both the auto and cross-corre- 
lation modes. The error calculation of Appendix A, however, shows that 
for small values of g the uncertainty involved in such a procedure makes 
the result almost meaningless; this situation evolves because the compu- 
tation requires the subtraction of two nearly equal numbers, each with 
non-trivial relative error. Although not used for computation the method 
1s outlined briefly below. 

With mirror surface, the received sound pressure is, as before, 
BP cos w(t-t'). Autocorrelation of this signal yields C(t): 


BP 20 sec 


C(t) = —-—=- i cos w(t-t') cos w(t-t'-t) dt 
20 0 


This function has a maximum at t = 0, equal to 


Care. 
B PA 20 sec 0 


Bene = eC = ——9°~ feeeeamnlnS maanlebestll) mackie ol 


I 
max 20 0 


d 


The surface is roughened and the autocorrelation repeated; the received 


Signal is ABP ,cos w(t-t') + p(t) 7 and its autocorrelation, evaluated 


inco 
at t=0 is D(t)|™”. 


aeBeP 20 sec 


D(O) = --+----- f cos¢ w(t-t") dt 
20 0 


cos w(t-t') dt 





70 0 incoh 





20 «0 f P incoh me 


The second term vanishes upon integration; the remaining two terms are: 


D(0) = AaB (1) : (rpcon) 
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Forming the ratio D(0)/C(0) gives a new quantity R: 


But B 


o> 
D(O) _ A'B Ld i Po 4 ne : aaa 
C(0) 
Bo / I a 
d d 


2 (14) = Citoieror) 


Cincoh ) 
e Ac —Nincoh 
Gail. 


Cn 
incoh a R - ac 


(Inirror) 


Determination of both A~ and R_ required the use of the correlator; 


2 


Since only one such instrument was available, the two quantities could 


not have been computed simultaneously, a fact which would have further 


added to total underta inty in the final result. 


D. DATA COLLECTION PROCEDURE 


A data run consisted of the following steps performed in sequence: 


1. 


With mirror surface and the transducers positioned at the desired 
grazing angle, the receiver and projector were physically aligned 
for maximum scattered return. 


Using pulsed sound, the ratio of unwanted direct-path pressure to 
mirror-surface-reflected-pressure was determined and recorded. 
Only data with direct path less than 3.0% of mirror reflection 
value were studied. 


With mirror surface and continuous-wave sound, the cross- 
correlation of driving voltage with received signal was measured; 
the value of the maximum of this function was read from an oscillo- 
scope after 180 seconds of correlator computation time 
(corresponding to six time constants of the integrators). 


The surface was then roughened, and the cross-correlation was 
repeated. The surface wave record, z(t), was tape recorded during 
this correlation, and the time of the correlation reading was 
marked by annotating the tape on a separate voice channel; duration 
of this record was three minutes. 


ae} 





5. At a later time the tape containing z(t) was replayed at original 
speed and autocorrelated with zero delay time. The resulting 
correlator output — mean squared wave height — was read and 
recorded when directed by the voice channel, corresponding to the 
time of previously having read sound correlation. 

The circuitry and equipment used in the foregoing steps are shown in 

figures 18 through 20. A separate data sheet was filled out for each 
data run, and the run was identified by its most important parameters; 
an example of run labelling is given in figure 21. 

Even though both RMS wave height and Na varied slightly with time 
throughout a data run (see section III-B-2-a) they were determined as 
nearly simultaneously as possible, each having been processed by inte- 
grators with the same time constant (the same integrators). It is believed 
that this method permitted the taking of meaningful data, the precision of 


which was not appreciably degraded by temporal variation of important 


parameters. 
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IV. ANALYSIS OF DATA 


A. DEPENDENCE OF COHERENT, SPECULAR SCATTER ON SURFACE 
ROUGHNESS FOR SIMPLE AND COMPOUND SEAS 


If surface shadowing is neglible, Beckmann and all others predict 





= a g : . ® 
Toon) 8 nan (Eq. 1). There is no reason, however, to believe 
that this function, derived for a Gaussian distribution of wave heights, 
would govern the scattering of specular, coherent sound from pure swell 
with its almost sinusoidal wave record (Appendix B). 


To compare eT 5 observations with theoretical predictions, 
I obs 
coh 


: - =e 
the ratio Io, (predict = 9 was formed for each data observa- 


tion. Departures of this ratio from unity can be attributed to one, or 
a combination, of the following factors: 


1. Non-validity of equation 1 in predicting the observed phenomenon. 





2. Experimental uncertainty (shown in Appendix A to be 27% for the 
worst case). 
3. Surface shadowing. 


Figures 22 through 25 show this quantity plotted against grazing 





angle for the various scattering surfaces. In all cases — with the ex- 
; . I op (obs) 
ception of figure 23 — the ratio Te predies). is within 8% of 
coh predic 


unity, well inside the worst case experimental uncertainty. Figure 24 
shows that coherent, Specular scattering from narrow-band swell alone is 


predicted by eS with at least the degree of precision found when the 


sea has a Gaussian distribution of wave heights. 





The curve representing 100 kHz and wind driven surface "C" shows a 
marked and significant departure from predicted behavior in the vicinity 
of 20° grazing angle. Three identical data runs were made on different 


days at 20° grazing to check the phenomerion's reproducibility (Figure 23); 
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all points were well above the predicted value of unity, and well beyond 
experimental error. This phenomenon was again manifest during the study 
of scatter for sound beam cross versus beam parallel to sea propagation 
direction (see figure 31). 

Since this gross departure from theory was noted only in the vicinity 
of 9 = 20° with the roughest wind-driven sea, these peyaneters were held 
constant, and the frequency was varied from 50 to 130 kHz in ten kHz 
steps. The resulting data points, plotted in figure 26, show the onset 
of the excess coherent component at just above 50 kHz. 


All data points for 20° grazina angle are a against surface 


I oni obs 
roughness, g, in figure 28. Here, I predict appears to be a 
Collie 


monotonically increasing function in the range 0 <gq< 3.0. The pAeno= 
menon cannot be attributed to increasing roughness alone, however, since 
at 6 = 30° even higher roughnesses were attained without this behavior, 
as shown in figure 27. Similar data, taken at 10° grazina anale, are 
shown in figure 29; the phenomenon in question does not appear, but it 
should be noted that for all points q< 0.4. Investigation of larger 
roughnesses at 10° grazing will require sound frequencies in the range 
130 kHz < f < 250 kHz, a region in which the EDO-327 transducers were not 
designed to operate. (250 kHz is the high frequency limit on Pues to 
the PAR-101 Correlation Function Computer). The study of this departure 
from theory is therefore incomplete, and investigation of larger rough- 


nesses (through higher frequencies) is necessary. 


B. COMPARISON OF CROSS-SEA WITH PARALLEL-SEA SCATTERING 
The width of the anechoic tank limited transducer separation to 52" 
(132 cm) for cross sea scatter (figure 4). If the distance from trans- 


ducer face to the beginning of the far field is taken to be x = 2a°/) : 


OF 


where a = piston radius = 1] cm, then x = 81 cm and 162 cm for 50 kHz 
and 100 kHz, respectively. Since the 52" spearation clearly placed the 
ensonified surface in the near field, for comparison it was considered: 
necessary to maintain this same geometry for parallel sea scatter (al- 
though much more space was available). 

The curves of figures 30 through 33, which compare cross and parallel 
scatter, differ = detail from the corresponding data taken at 18' trans- 
ducer spacing, possibly due to operation in the near field. For this 
reason, the data obtained at 52" separation were not used to study depend- 
ence of specular, coherent scattering on roughness, but were only used to 
infer information concerning cross versus parallel] scatterina. 


] tops) 
CO ‘ 
Figures 30 through 33 show ee predict) versus arazing angle 


for both cross and parallel scatter. In almost every case, it is seen 
that cross scatter gave a somewhat stronger specular coherent return than 
was researc for the corresponding parallel case. The noticeable drop 
which occurred at 100 kHz, 30° grazing angle was outside the bounds of 
experimental uncertainty; since it was not present at the greater trans- 
ducer separation, it appears to be attributable to either near field 
operation or to the reduction in ensonified area (by a factor of about 6). 

Tt should be noted that the larger than expected coherent component 
at 0 = 20°, sea "C" occurs. in figure 31, for both cross and parallel 


scattering. 


C. SURFACE SHADOWING 

The shadowing function, S(6), derived by Wagner [Ref. 5] and defined 
in section II-B uses ray acoustics, neglecting diffraction as well as 
multiple reflections. The ray approximation should improve with increasing 


o/X , becoming precise in the limit as o/d>~ . Wagner states that the 
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correct use of the shadowing function for specular, coherent intensity is 


multiplication by s“(e), Thus 


I | (obs) 
oo = S*(0) & ip Sy T_ (predict) *(8) 
mirror Loh 
/ 
The arguments of S(6) are 6, 6, and L where = surface correlation 


distance in the scattering direction. The sound wavelength, », does not 
appear, but use of the ray theory implies that o >> A. 

The spatial surface correlation function measured in the scattering 
direction (up-down wind) resembles a damped cosine curve [Ref. 7]. If 
‘the correlation distance (distance for the correlation function to reach 
Ve of its initial value) is determined using the decay envelope, L_ 14.0 
cm; however, if L is computed from the cosine-like function in its first 
quarter period (prior to the first crossing of the zero correlation axis) 
L 3.0 cm results. The correlation distances are indicated to be approxi- 
mate, since they were obtained for a slightly smaller sea generated by a 
set of fans wiffch is no longer iristalled. 

The transducers used limited the value of o/a to 0.3 and less. 
Although this obviously does not satisfy the requirement that o >> id, 
data obtained were coared to Wagner's predictions (Fig. 34) plotted for 
convenience at o/i = 10. Lack of experimental data in the region 0.5 
0.5 << ofa < 10.0 precludes extrapolation and leaves unanswered the 
question whether the correlation distance computed from the decay envelope 
(L = 14.0 cm) more accurately predicts the degree of shadowing than does 


= 3.0 cm 
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V. CONCLUSIONS 


The unpredicted growth of the coherent component of scattered inten- 
sity in figures 23, 26, and 28 may be the same phenomenon described by 
Medwin [Ref. 1]. While conducting normal-incidence, specular scatterina 
from a wind driven surface, Medwin observed that above a critical fre- 
quency "the instantaneous scattered intensity is a nonstationary, oscil- 
lating function of frequency, varying about the asymptotic value." He 
postulated that this oscillation with frequency was due to alternate 
destructive and constructive interference of incident and reflected sound 
waves from Staeteui cat prominent facets of the surface, whose distri- 
bution of wave heights was Gaussian. Medwin suggested that the dimension 
causing the interference was perhaps the distance between mean positive 
and mean negative wave heights, and noted that the phenomenon was absent 
for o << 4. He defined the critical frequency as that for which: 1/4 = 
facet peeranion (the first interference with increasing frequency). 

Non-normal sound beam incidence and surface anisotropy greatly compli- 
cate definition of the facets responsible for interference. Examination 
of figures 30 through 33 (cantnan Tardy the 100 kHz curves of Fia. 31) 
indicates little, if any dependence of the behavior on surface anisotropy; 
thus, if this is an interference phenomenon, it suggests that the. qover- 
ning parameters are functions of 4, o, and 9, and are not influenced by 
correlation distance. | 

Comparison of parallel-sea with cross-sea scattering shows that the 
cross-sea geometry gives a coherent component which is in general consis- 
tently larger than that obtained for the corresponding parallel-sea case 


(Figs. 30 through 33). This difference in coherent component may have 
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been due to the different off-specular rays for the two cases; these 
rays, which are excluded in theoretical solutions, are inevitably collected 
during an experiment in spite of the use of directional transducers. 

Specular, coherent scattering from pure swell (whose wave height 
probability density is strongly non-Gaussian) appears to closely obey 
the e 9 prediction, although this function was derived for a Gaussian 
distribution of surface wave heights. The high-frequency limit of the 
ED0-327 transducers and the relatively small RMS swell wave height 
(o < 0.15 cm) restricted the attainable values of g to 0<g<0O.1. 

For completeness, behavior of specular, coherent scattering in the range 
0.1 < g < 5.0 should be studied. 

The energy spectral density for sea-swell combination "E" [Fig. 15] 
shows that about 20% of total energy is carried by the swell. The corre- 
sponding wave height probability density function, however, appears as 
near to Gaussian as any of the purely wind-driven seas (Fig. B-5). Non- 
linear wave interaction caused the energy spectrum of the sea-swel | 
combination to be significantly different from that produced by super- 
position of its two constituent spectra (see Figs. 13, 14, and 15). 

Wagner's Shadowing Correction (Sec. II-B), based on ray acoustics, 
tended to overcorrect, due no doubt to the partial illumination of 
"Shadowed" areas caused by wave diffraction. In figure 34 it is not 
possible to extrapolate to a higher value of o/d, where the ray theory 
should be more valid, in order to determine whether the proper correlation 
distance to use in the computation of Wagner's correction should be 
derived from the envelope of the spatial correlation function. 

At the time of this experimental work, transducers were being fabri- 
cated to cover the frequency range 120 < f < 240 kHz. Investigation of 
specular scattering in this frequency range is necessary to clarify many 


of the phenomena observed. 66 





APPENDIX A — ANALYSIS OF ERROR 


Several important calculated quantities concerning scattered sound 
depend on oscilloscope readings and the wave staff calibration. It was 
possible to estimate the degrees of uncertainty in these constituent 
quantities, and their effect on the relative coherent intensity and 
incoherent intensity compared with theory is now considered. The method 
used is the finite differential approximation to the derivative. 

| | | /1 , obs) 
. Uncertainty in \ Toon PredTet] 
This ratio, which compares the observed to predicted specularly 


scattered coherent intensity, is called I for brevity in the analysis. 


oo I on obs) : wie ; 1209 

= T_.. (predict) = = = 
oh predict ag 

_ NE ok%0%(2 sin @)* 
Me 

where: N = maximum value of rough surface cross-correlation, 
read from scope 

M = maximum value of mirror surface cross-correlation, 


read from scope. 


It is assumed that k and 6 are known exactly; sources of uncertainty are 


M, N, and o. 

Let: sly 
ss = relative uncertainty in I due to N 
bly 
Te = relative uncertainty in I due to M 


6/ 


=a = relative uncertainty in I due to o 


The total relative error is the sum of the contributions: 


cleotay Ne. SNe ol 


I I 


By partial differentiation: 
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(3) 3L = pe a? 5 (2 sing)? eF = 29L 
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Using the finite differential approximation and rearranging, equations 


(1), (2), and (3) become: 
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where ae ‘ a ,» and ~ are the constituent relative uncertainties. 


It can be seen taht the "worst case" will occur when the following three 


conditions are simultaneously satisfied: 
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a. N is small 
b. gis large 


c. 6M is negative, while 6N and 6&o are positive. 


It should be noted that the value of M is controllable; its magnitude is 
set for almost full oscilloscope deflection at the start of each data run 
by adjusting the correlator input gains. 

For the "worst case" analysis, typical parameters from a large rouah- 
ness data run were used. In all cases, the oscilloscope readinas were 


judged accurate to f.1 cm. 
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If the method of section III-C-2 had been used to isolate the 
incoherent intensity, uncertainties of unacceptable magnitude would have 


resulted, as will no be shown. Define: 


me 


yg ue, ee (from section III-C-2) 
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where: D = maximum value of rough surface autocorrelation, 
read from oscilloscope 

C = maximum value of mirror surface autocorrelation, 
read from oscilloscope 


By the method of finite differentials 


6D 6C 2 6M SN 
Hoel . 2) ia oe ee 
a om 


R-A 


Using parameters which occurred during an actual data run, (Run 25-E-100): 








G = 56 
R = ./713 
Ao = 671 
6C _ 6M _ - .1 cm a é | 
Se a en (negative for worst case) 
oo oN os 1 cm e 
_— oe a 


Or 
| a 


- (.713)(.03) + 2(.671)(.03) 2 4.45 = 145e 


The difficulty arises in the term (R - Ac); the two numbers are nearly 
equal, and each has significant relative error. The above calculation 
gives a conservative estimate of the uncertainty, since it does not con- 
sider the non-simultaneity of cross and autocorrelations in the presence 
of statistical temporal variations. In order to obtain the incoherent 
intensity, it would have had to be calculated or measured directly, rather 


than computed by subtraction. 
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APPENDIX B — PROBABILITY DENSITY FUNCTIONS 
OF THE SURFACE WAVE HEIGHTS 


The probability density of wave heights was determined for each of 
the representative scattering surfaces as follows: 

1. With the wave staff located at the center of the scattering area, 

a ten minute tape recording was made of the wave record, z(t). 

2. The wave record was replayed at ten times speed into the analog- 
to- digital converter of the CI5000/SDS9300 hy brid computer. The converter 
sampled the analog signal, z(10t*), at a rate of 333 Hz, until an ensemble 
of 10" digital wave heights was generated. (t* is time in the playback- 
digitize frame). 

3. The mean and standard deviation of the ensemble were computed; the 
ensemble mean was then adjusted to zero by subtracting the previously 
calculated mean from each ensemble member. 

4. The SDS9300 digital computer then sorted the ensemble values into’ 
intervals of width co, where o = ensemble standard deviation, and c_ was 
a constant to be optimized. 

5. Total "area" under the resulting probability density histogram 
was 10° (the number of observations). Dividing each ordinate value by 
10° then gave a probability density whose area was normalized to unity. 

The appearance (and utility) of the density functions obtained was 
dependent upon the interval width selected; larger widths produced a 
smooth plot, but tended to suppress all but the gross features, while too 
narrow an interval gave an irregular plot due to sparsity of interval 
population. Figure B-2, made with an interval width = 0.050 , illustrates 


the latter phenomenon. 
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For all wave-height probability densities (except pure swell), a 
standard Gaussian curve normalized in units of o was plotted for compar- 
ison. These analytical curves are represented by solid lines, while 
the observed density function is given by points. 

The wave-height probability density of pure swell (Fig. B-4) 
resembles that of a sine wave, even though its frequency-domain 


structure (Fig. 14) is by no means a delta function. 
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